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ABSTRACT 

A new method is presented of determining the three-dimensional sound-speed structure and flow velocities in 
the solar convection zone by inversion of the acoustic travel-time data recently obtained by Duvall and coworkers. 
The initial inversion results reveal large-scale subsurface structures and flows related to the active regions, and 
are important for understanding the physics of solar activity and large-scale convection. The results provide 
evidence of a zonal structure below the surface in the low-latitude area of the magnetic activity. Strong converging 
downflows, up to 1.2 km s ', and a substantial excess of the sound speed are found beneath growing active 
regions. In a decaying active region, there is evidence for the lower than average sound speed and for upwelling 
of plasma. 

Subject headings: convection — methods: data analysis — Sun: activity — Sun: oscillations — sunspots 


1. INTRODUCTION 

Sound waves excited in the subsurface granulation layer 
propagate into the Sun’s interior along ray paths which are 
relraeted toward the surface because the speed of sound 
increases with the depth. When the sound speed becomes 
sufficiently high the waves are reflected and return to the 
surface at some distance from their origin. The travel time of 
the waves depends primarily on the speed of sound (wave 
group velocity, more precisely) and on the velocity of flow 
along the ray paths. The effects of the sound-speed structure 
and of flows can be separated by measuring the travel time of 
waves propagating in opposite directions along the same ray 
paths. Duvall et al. (1993) have discovered that the travel time 
can be measured for groups of rays despite the stochastic 
nature of solar sound and turbulent noise on the surface. Such 
measurements provide a means for probing the three-dimen- 
sional structures of the internal sound speed and flows. 

Recently, Duvall et al. (1996) have published maps of the 
travel times of sound waves propagating in local regions in the 
convection zone. The maps represent measurements of the 
time tor acoustic waves to travel between points on the solar 
surface and surrounding annuli. The measurements were 
obtained by estimating cross-correlation functions of fluctua- 
tions of brightness of the solar surface, caused by 5 minute 
acoustic oscillations. The waves are stochastically excited by 
granulation uniformly distributed on the surface; therefore, 
there are waves propagating between any two points on the 
surface. However, because of the high level of noise, the data 
are averaged over large surface areas. Typically, the measure- 
ments represent a mean travel time between a central point 
and surrounding annuli, which were chosen in four radial 
distance ranges, A: 2.5-4.25, 4.5-7, 7.25-10, and 10.25-15 
heliographic degrees (1° « 12.15 Mm). The acoustic rays that 
propagate to these radial distances span the upper 30% of the 
convection zone, penetrating to a depth, z, of about 64 Mm 
(Fig. 1). 

Two sets of data were obtained during 17 hr of continuous 
observations at the South Pole in 1991. The first represents the 
mean of the travel time r for the rays propagating from a 
central point to the surrounding annuli and the time r for the 
rays propagating in the opposite direction. The second set is 


the difference of these travel times. The mean travel time, 
0.5(7* + t ), depends primarily on the sound speed along the 
ray paths, while the difference 7 4 — t is sensitive to the flow 
velocity. The observed area is shown in Figure 2 (Plate L8). 
The spatial resolution of the data is 1 ? 75 (21.3 Mm) in both 
longitude and latitude. 

We have developed an inversion method to infer the 
structure of the sound speed and flow velocities from the data. 
The technique is somewhat similar to seismic tomography 
(c.g., Howe, Worcester, & Spindel 1987; Iyer & Hirahara 
1993); however, it is more complicated computationally be- 
cause each data point is no longer a one-dimensional line 
integral but is an average over ~1(T ray paths that sample a 
large three-dimensional region (Fig. 1). 

2. THE INVERSION METHOD 

The travel time, r n along a single ray path, F,, in the 
geometrical optics approximation can be expressed in terms of 
the sound speed (wave group velocity), c(x, /), and the flow 
velocity r(x, /): 


Ti(t) 


i. 


ds 

c(or, t) + rU, /) • n 


( 1 ) 


where x is the spatial coordinate, / is the time, .v is the distance 
along the ray, and n is a unit vector tangent to the ray. The sign 
of v • n depends on the direction of propagation; therefore, the 
travel time in opposite directions differs due to the effects of 
flows. Since variations of the travel time in the area of 
observation do not exceed 5%, the effects of flow' and sound- 
speed inhomogeneity can be considered, in the first approxi- 
mation, as perturbations to a horizontally uniform static 
reference model with a sound speed c n ( jr). Then the reference 
travel time is 7 (l , = ^ dslc u (x). 

Fermat’s principle applied to the variation of the travel time 
gives 


8t,' (/) = 7 , + (/) - 7 n , 


8c(j r, /) ± r(x, /) • n 
ci(x) 


ds , 
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p t(i ] Regions of ray propagation (shaded areas) corresponding to the 

measurements by Duvall et al. (19%), as a function of depth, z, and the radial 
distance, A, from a point on the surface. The rays are also averaged over 
circular regions on the surface, forming three-dimensional figures of revolu- 
tion. The dashed lines are the boundaries of the cells of the inversion model. 

where the plus sign is for the rays propagating from the central 
point, A = 0, to the annuli, and the minus sign is for the rays 
propagating in the opposite direction. The average and the 
difference between the travel time residuals are 

1 f 8c(x,t) . 

8rr(t) = 2 [ 6t '* (,) + 5t . ] = - ds ' ( 3 ) 

* I’ll./ 

r(x, t)*n 

t) = St/ (t) ~ St, = -2 ( .. (jr p ds. (4) 
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We assume that 5c and v did not change during the observa- 
tions and represent them by a discrete model. In the model, 
the three-dimensional region of wave propagation is divided 
into rectangular blocks. We also assume that the relative 
perturbations of the sound speed, 5 c/c, and the ratio of the 
flow velocity to the sound speed, M = vie, are constant within 
the blocks. Then 

5 In c(x, t) = 5 In c ijk , M(x,t) = M ljk , (5) 

where 5 c ijk and M lJk are the unknown parameters of a block of 
the spatial coordinate indices i,j, and k. 

According to Duvall et al.’s (1996) averaging procedure, the 
travel time measured at a point on the surface is the result of 
the cumulative effects of the perturbations in each of the 
traversed rays of the three-dimensional ray systems (Fig. 1). 
Therefore, we average equations (3) and (4) over the ray 
systems corresponding to the four radial distance intervals of 
the data, using approximately the same number of ray paths as 
in the observational procedure. As a result, we obtain two 
systems of linear equations that relate the data to the sound- 
speed variation and to the flow velocity: 


st) v ;; = X A T 8 ln c,A i 

ijk 

(6) 

8rt‘.= X 

( 7 ) 


ijk . « 


where A. and \x label the central points (A = 0) of the ray 
systems in the observed area; v (1, ■ . ■ , 4) labels the radial 
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Pi (; — Azimuthally averaged perturbation of the sound speed as a 
function of latitude in the upper layer, 0-16 Mm (solid curve), and in the second 
layer, 16-32 Mm (dotted cun'e). 

distance ranges (or annuli); /, j, k are the indices of the 
rectangular blocks of the three-dimensional structure model; 
and a denotes three components of the flow velocity. Equa- 
tions (6) and (7) have been solved by a regularized least-squares 
technique (Paige & Saunders 1982). The mathematical details of 
the inversion procedure and tests will be published elsewhere. 

3. THE INVERSION RESULTS 

The travel-time anomalies have been measured by Duvall et 
al. (1996) for four sets of rays propagating from each of 
72 X 52 points to the surrounding annuli (Fig. 1). The points 
are equally spaced in the observed area of 123° X 106 on the 
solar surface (Fig. 2). For inversion of the data, we have 
adopted a four-layer discrete model. The layers which all are 
16 Mm thick have been divided into 72 X 52 rectangular 
blocks of approximately the same size as the spacing between 
the data points. 

The results of inversion are shown in Figure 2 and in Figure 
3 (Plate L9) as maps of the sound-speed perturbation, Sc/c, 
and the three components of the flow velocity, v, in the four 
layers of our discrete model. They can be compared with the 
Ca + K-line intensity map in Figure 2. Perhaps the most 
interesting detail on the sound-speed map is the equatorial 
belt of the lower sound speed in the upper layer and higher 
speed in the lower layers. The belt is surrounded by narrow 
regions of the higher sound speed at latitude «±15°. This 
could be evidence of a thermal shadow due to a band of 
magnetic field in the convection zone. The thermal shadow 
was suggested by Parker (1987), who argued that the broad 
bands of azimuthal field in the convection zone should pro- 
duce cool shadows which press the field downward, and force 
it to remain in the lower convection zone, balancing the 
magnetic buoyancy. Evidence for the thermal shadow has been 
previously reported by Kuhn, Fibbrecht, & Dicke (1988) from 
photometry observations of the surface brightness, and by 
Goode & Kuhn (1990) from splitting of p-mode frequencies. 
Figure 4 shows latitudinal variations of the azimuthally aver- 
aged sound-speed structure in the upper two layers. In the 
upper laye r, the averaged relative variation of the sound 
speed, Sc/c, is ~-5 X 10 \ the magnitude of which is some- 
what higher than the variations of the surface temperature 
(Kuhn 1988), suggesting that the main changes occur beneath 
the photosphere, as was also noted by K uhn ( 1988) and by 
Goode & Kuhn ( 1990). In the second layer, Sc/c ~ +6 X 10 5 . 
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The variations of the sound speed can he due to temperature 
change or due to magnetic field: 

_ 8T B 1 

c 2 T H irp ’ 

where p is the gas pressure. If the inferred variations are 
entirely due to magnetic field, then B ~ 5 x 1() 4 G for the 
mean pressure in the second layer, P ~ 10" dyn cm 2 c ! . This 
field strength is well above the low limit of 3 X 1() 3 G 
established by Parker (1987). 

It is interesting to note that the boundaries of the belt 
(«±10°) coincide with the most active latitude at the time of 
the observations (Solar-Geophysical Data 1991). There is also 
a hint for two narrow strips of reduced sound speed at latitude 
±22°, where active regions also emerged frequently. 

The perturbations of the active regions can be seen in the 
two upper layers, but not deeper. The velocity maps (Fig. 3) 
show strong downflows of « 1 km s 1 around the sunspots, and 
some much weaker isolated upflows at the boundaries of the 
equatorial belt. The spatial resolution of the data is insufficient 
to study the structure and dynamics of individual sunspots. 
However, it is interesting to note that the peak of the 
downdraft velocity is observed between the spots in the bipolar 
group NOAA 6432 located at latitude - 10° and longitude 20°. 
This group results in two separate strong negative sound-speed 
perturbations in the top layer (0-16 Mm), which merge into a 
single perturbation in the second layer (16-32 Mm). Figure 5 
(Plate L10) (top) shows a vertical cut through this active region 
at a fixed latitude. The color map represents the sound-speed 
structure, and the arrows the projection of the velocity field 
into the longitude-depth plane. The convergent character of 
the flow is evident. This active region was only about 5 days 
old. The other young active regions in the observed area seem 
to show a similar pattern. In contrast, the much older decaying 
region NOAA 6427 (Fig. 5, bottom) shows reduced sound 
speed and diverging upflow. Another active region, NOAA 
6437, with similar characteristics disappeared a day after it was 
observed. The flow patterns in the active regions seem to be 
consistent with the Doppler observations of the plasma motion 
on the solar surface (Gopasyuk 1967; Zwaan 1985). 

In the flow maps (Fig. 3), the largest flow velocities are 
observed at the active latitudes ±10° and ±22°. Some areas 
show strong flows without any apparent activity. At least, some 
of them could be the places where new magnetic structures 
emerged later— e.g., the area of strong converging flow seen in 
the two upper layers at latitude -8° and longitude 10° was 
close to the place where a new large group of sunspots, NOAA 
6440, was born 2 days after. 

4. DISC USSION 

Time-distance helioseismology pioneered by Duvall et al. 
(1993) gives us a three-dimensional look into the Sun’s 


interior, providing unique information about the structure and 
flows associated with magnetic field and turbulent convection. 
The interlink between convection and magnetic field results in 
the remarkable phenomenon of the solar cycle which has no 
analogy in the terrestrial laboratory. 

In this Letter we have presented a technique and some 
initial results of inversion of the travel-time data which are 
sensitive to the sound-speed structure and to the flow velocity 
along sound-wave ray paths. A three-dimensional inversion 
method based on Fermat’s principle and a regularized least- 
squares technique have been applied to infer the sound speed 
and the velocity of flows from the observations obtained at the 
South Pole on 1991 January 4-5 by Duvall et al. (1996). The 
spatial resolution of the inversion is 1 ?75 in both longitude and 
latitude and 16 Mm in depth. The results for a four-layer 
model reveal large-scale subsurface sound-speed structures 
and flows related to the active regions. 

We have found evidence for a large-scale low-latitude 
azimuthal structure which, possibly, can be associated with 
azimuthal magnetic flux bundles and a thermal shadow above 
them (Parker 1987) and, also, with thermal banding of solar 
convection (Brummell, Cattaneo, & Toomre, 1995). The mag- 
netic flux bundles probably give rise to the magnetic U-loops 
emerging through the surface and observed as groups of 
sunspots (Parker 1994). However, the maximum amplitude of 
the azimuthal perturbation is found in the second layer, from 
16 to 32 Mm, whereas, from theoretical arguments, the flux 
bundles should be concentrated near the bottom of the 
convection zone. The time-distance seismology has not yet 
reached that depth. Also, the current evidence is rather weak 
because of the large observational errors, but future travel- 
time data will certainly provide the information about the 
large-scale azimuthal structure, important for testing dynamo 
theories. 

The emerging H-loops are expected to result in upwelling of 
the plasma in the active region. Instead, we see converging 
downflows in the places of recently emerged groups of sunspots 
and upwelling in decaying active regions, in agreement with the 
flow pattern observed on the solar surface. Also, there is evidence 
of different sound-speed structure in the developing and decaying 
regions. Obviously, continuous observations are necessary to 
study the dynamics of the active regions in detail. 

Long-term tomographic monitoring of the solar convection 
zone using high-resolution data from the MDI experiment on 
SOHO (Scherrer et al. 1996) will provide a means for under- 
standing the physics of solar convection and activity. 

1 thank Tom Duvall, Gene Parker, and Juri Toomre for 
useful discussions. This work was supported by a grant from 
NASA. 
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